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Abstract 

The aim of this paper is twofold. First, we obtain the explicit exact 
formal solutions of differential equations of different types in the form with 
Dyson chronological operator exponents. This allows us to deal directly 
with the solutions to the equations rather than the equations themselves. 
Second, we consider in detail the algebraic properties of chronological 
operators, yielding an extensive family of operator identities. The main 
advantage of the approach is to handle the formal solutions at least as 
well as ordinary functions. We examine from a general standpoint lin- 
ear and non-linear ODEs of any order, systems of ODEs, linear operator 
ODEs, linear PDEs and systems of linear PDEs for one unknown function. 
The methods and techniques involved are demonstrated on examples from 
important differential equations of mathematical physics. 

1 Introduction 

There are several universal approaches to solving DE. The most promising of 
them are Lie symmetry methods m, |2] and advanced methods for the linear 
equations which are based on differential Galois theory 

The aim of this contribution is to show examples of the possibilities for solv- 
ing differential equations with the method based on Dyson's operator solution 
representation in the form of time-ordered exponentials |1] (see also It 
is well-known that Dyson's use of time-ordering is the fundamental conceptual 
tool in quantum field theory. This tool has now become a natural part of many 
branches of physics and is even used in parts of engineering. 

The method allows us to find the explicit exact operator solutions of many 
problems if we succeed in reformulating a given problem into a first-order linear 
one. The next step of the approach is directed to transforming the operator 
solution into a more practical, calculable or useful expression. The use of oper- 
ator methods is efficient in manipulating time-ordered expressions with help of 
exponential identities like the Baker-Campbell-Hausdorff formulae [Z|, [H|. Our 
second aim in this paper is to give general methods for obtaining such identities. 
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Wc have in mind that at every stage the chronological exponential can be ex- 
panded in formal operator power series neglecting all the convergence problems 
that can appear with, for example, analytical functions. 

As a result, this method allows us to deal directly with the solutions to the 
equations rather than the equations themselves. 

The main advantages of the method are in its compactness, clarity and 
simplicity. It is also essential that wc can handle the operator expressions in 
manner similar to ordinary functions. 

It is clear that there are some correlations between the operator method and 
other common methods, especially the Lie approach. 

We choose example equations for the operator method from important dif- 
ferential equations from many branches of science. The main aspects of the 
method arc quite easily generalized for more complicated equations. 

As the paper is addressed to practical problems in many branches of physics 
and engineering for which solutions are not yet at hand and to avoid complica- 
tions, we choose a simple (and non-rigorous in places) manner of exposition in 
that we can always verify the solutions by substitution into the initial equations. 
Analytical restrictions are obvious from the context as a rule. 

2 Definitions and notations 

An operator A is defined as a mapping of a function / (from a ring R) into a 
function A/ (from the same ring R). An operator A is linear if it maps any 
two functions (j) and ip in such a way that 

A{(p + = A(f) + Atjj. 

Operators will be written in the bold font seen above to avoid any ambiguities 
with the exception being the usual differential operator notation. In every 
case throughout this paper, we shall be dealing exclusively with operators that 
are linear. 

The operator A is derivative if 

A(A + B)/ = AA/ + AB/, 

AAf = (AA)/ + AAf 

for any / (from a ring R). It follows from this definition that the inner derivative 
of an operator A is 

(AA) = AA - AA = [A, A], 

where [A, A] is the conventional notation of the two operator commutator. In 
this paper we will consider only differentiable operators A in sense that the 
commutator [A, A] exists. 
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If A = ^ then such a definition of the inner derivative of the differentiable 
operator A.(t) is in accord with the classical analytic definition 

dA{x) _ ^.^ A{x + e) - A(a;) 

dx e^o e 

We will find need for an operator-valued function F{A) which is an operator 
itself. Only functions which can be described by power series are considered in 
this work. 

Important examples for this paper are the following types of linear operators; 

where /, are arbitrary functions of t,xi, ...,Xn, Li is obviously a derivative and 

d d 
L2 = f(t,a;i,...,a;„,— — ). 

Though L2 is a generalization of Li, nevertheless L2 is not a derivative in 
general, as e.g., {-§^Y = does not satisfy the above definition of derivative. 

3 The linear first-order differential equations and 
basic identities for chronological operator ex- 
ponents 

3.1 Dyson's form of the solution of first-order linear dif- 
ferential equation 

Let us start with the first-order linear differential equation with respect to t for 
u{t,p) = u{t) : 

^=L(*M*) (2) 

with the initial condition 

u{t)\t=a = V, (3) 

where p is a set of parameters say xi, Xn, L(i) = hit, p) is a linear operator 
which does not depend on ^ explicitly, and v is some function of p. 
This equation can be solved by the following iteration scheme 

UQ{t) = V, 

9u„{t) 

—Qf - = Ht)Un-l{t), 

which for t > a leads to 

u{t) ='E,{t,a;L{t))v, (4) 
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where operator E has the foUowing series expansion (1 is the identity operator) 
E(t,a;L(t)) = 1 + / dTL(T) + 

J a 

■■■+ dn dT2... dTn L(ti) L{t2) . . . L(r„) + . . . (5) 

J a J a J a 

The solution Q (taking into account jSJ) has the fohowing properties: 

(a) . Operator E satisfies the operator differential equation 

^ = Mm (6) 

with initial condition 

nt=a = 1, (7) 

or the equivalent integral equation 

E(i) = l+ / dTL(r)E(r). 

J a 

(b) . All operators L in the series expansion of E are ordered in sense that 
any operator term L(ri) L(t2) . . . L(r„) corresponds to the requirement that 
Ti > T2 > ■ ■ ■ > Tn- The ordering here is very essential as operators L(Ti) and 
L(tj) do not commute in general. 

The property (a) suggests that the operator E is a kind of exp{/^* dr L(r)} 
but expansion of the exponent here does not lead to the required operator 
ordering (b), as for example, the third term of the E expansion 



f dn f ' dT2 L(Ti)L(r2) 

J a J a 



differs from the third term of the exponent expansion 



i^*dTi£dT2L(Ti)L(T2). 



The next trick was proposed by Freeman J. Dyson [1] (see also 0, El)- Let 
us demonstrate it on the third term of the E expansion. By changing notation 
of integration variables and integration order we obtain 



/ / 'dT2L(ri)L(T2) = j dn f dT2L(r2)L 



so 



f dn /" 'dr2L(Ti)L(T2) = i / dn{ r dT2Uri)Mr2)+ I dra L(t2) L(ti)}. 
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If we now introduce a T-operator defined as 

T{L(ri)L(r2))} = 
then n-th term of the right-hand side of ISJ) can be transformed to 



L(ti)L(t2)) if ri > r2; 
L(t2)L(ti)) if ri < r2, 



/ dn dT2... dTn L(ti) L(r2) . . .L(r„) 

./a ./a ./a 

-fdnfdT2...f dTn T{L(ti) L(r2)...L(r„)}. 

J a J a J a 

Removing the T-operator outside of the integral signs means we can express 
E(t,a;L(i)) = T{1+ /" dTL{T) + 

J a 

■■■ + ^ dTiJ^ '^^^ ■ ■ ■ / "^^^ ^^^^^ ^^^^^ ■ ■ ■ ^ ■ ■ ■ ^ 

or in the final form with the chronological operator exponential 

E{t,a;L{t)) = Texp{ f dTL(T)}, (8) 

J a 

where we have in mind that IHI represents the series expansion and that 
T{L(ri) L(t2)...L(t„)} =L(r„J L(r„ J . . . L(r„ J. 

The expression © is the formal solution of the operator differential equation 
(O with initial condition (TJ. The formal solution of the problem @ is 
obviously 

u{t)^Texp{f dTL{T)}v. (9) 

J a 

The main advantage of solutions in the chronological exponential form is the 
explicit dependence of all parameters of the problem and the relatively clear 
way to obtaining approximate non-formal solutions through series expansion 

Texp{/" dTL{T)} = l+f drL(T) + 

J a J a 

■■■+ dn dT2... dTn L(ri) L(t2) . . . L(T„)-h 

J a J a J a 

f dTi r dT2... r dTn+1 L(ri) L(r2) . . . L(7:„+i) Texp{ ( d^UO) ■ 

J a J a J a J a 
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Analogously, the solution of the problem Q for t < a is similar to 0, 
but instead of the operator T it is necessary to replace it with the operator Tq , 
which represents a product of L's, in reverse order 

To{L(ti) L(t2)...L(t„)} = L(t„J L(T„J...L(r„J. 

The problem ||2Jl, © is deterministic in the sense that if we start from 
an initial condition at t ^ a we can find a unique u{t, p) for any t. If we 
now consider equation (|2Jl with a new initial condition u{t^p)\t=b = "Viip) and 
recalculate M(t, p) for t > b we obtain (for a < b < t) 

ft pt ph 

u(i) = Texp{/ dTL(r)}t; = Texp{ / dT'L{T)}T ex:p{ l drL(r)}?; 

J a J h J a 

for any v. The operators T here and everywhere in this paper act on operators 
of their exponentials only. That is, for a < b < t we have the following operator 
identity 



pt pt pb 

Texp{/ dTL(r)} =Texp{ / dr L(r)} T exp{ / rfTL(r)}. 

-/a 'J b J a 



(10) 



Analogously, if on the second stage we recalculate u(t,p} from f to a back- 
ward, we obtain that 

pa pt pt 

V = Toexp{ / dT'L{T)}u{t) = Toexp{- / dTL(r)}Texp{ / dTL(r)}z; 

J t J a J a 

for any v. Therefore 

Toexp{- /" dTL(r)}Texp{ /" dTL(r)} = l. (11) 

J a J a 



So the operator 



is inverse to 



Toexp{- f dTL(T)} 

J a 

Texp{/" rfTL(T)} 

J a 



and these operators commute. 

It can be easily seen that the operator 



E-i(t) -Toexp{ / drUr)} (12) 

J a 



is the solution of the following operator differential equation 



=E-i(t)L(t) (13) 
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with an initial condition 

E(t)|t=a = 1 (14) 

and that ^ 

E-\t) = l+ f dTE-\T)L{T) (15) 

J a 



which leads to the following series expansion 
ft 



E-i(t) = Toexp{ /" dTL(T)} = 1 + /" dTL{T)+ 

J a J a 

\- dTi dT2... dTn L(r„) L(t„_i) . . .L(ri)+ 

'J a 'J a <J a 

I dn f ' dT2... f " dT„+iToexp{ f deL(e)}L(T„+i)L(r„) . . .L(ti) . 

J a J a J a J a 

(16) 

It follows from properties H13|) - (|16|l (or from (|10|l ) that fov a < b < t 

pt nb pt 

To exp{ / dT L(r)} = Tq exp{ / dr L(t)} Tq exp{ / dr L(r)}. (17) 

J a J a J b 

With property © at hand one can easily verify that the solution of the 
inhomogeneous linear differential equation 

du(t) 

— ^ = L{t)u{t) + Cf>{t) {u{t)\t=a = V) 

at 

has the following operator form [t > a) 

M(t)=:Texp{/" dTL{T)}v+ f dr Texp{/ d^ L(e)}<?!)(T). (18) 

J a J a J r 

The operator technique in the form given above is immediately suitable only 
for solving the linear problems. As it will be demonstrated later, an analogous 
approach is fruitful for more complicated problems. However, before touch- 
ing such problems we need to consider the algebra of chronological operator 
exponentials in more details. 

3.2 Basic transformation identities for chronological op- 
erator exponents 

For any linear operator L(t) (not necessarily invertible) we can form the chrono- 
logical operator exponential {here and later on t > a) 

Texp{/" dTL(T)}, 

J a 

which is always invertible and differentiable with respect to t. 
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Suppose we have an arbitrary linear invcrtible difFerentiable operator A{t). 
From the obvious identity 

^ = {^A-\t)}A{t) (A(0|.=. = A(a)) 

we obtain ^ 

A(t) = Texp{^ dr^^A-i(T)}A(a). (19) 

It follows that the chronological operator exponential can represent the t-dependence 
of any linear invertible differentiable operator, and each such operator corre- 
sponds to the following linear first-order differential equation 

^ = {^A-(t)Ht) ium..^.) 

with the solution 

u(t) = A{t)A-\a)v. 

As the product of any two of such operators is the linear invertible differen- 
tiable operator (i.e., they form group) then the set of all chronological operators 
forms group too. 

If we invert (|19|) 

A-\t) = A-i(a)To exp{- £ dr A-i(r)} 

and replace A^^{t) by B(t) we get 

B(t)=B(a)Toexp{-y dr ' B(r)} ^ 

B(a)Toexp{ f drB-^r)^^^} 

and conclude that any linear invertible differentiable operator can be represented 
as an inverse ordered operator exponential too. Hence, the two main forms of 
opposite ordered operator exponents can be expressed through each other (see 
iPI). below). 

Let us consider the following operator 



K(t) = b(t)Texp{ /" drA(r)}, 

J a 



where b(i) is a linear invertible differentiable operator and A{t) is a linear one. 
Differentiating K with respect to t we have 

9K|i) ^ ^9b^b-i(0 +b(OA(Ob-i(t)}K(0 (20) 
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and 

K(i)|t=a =b(a). (21) 

As it was shown in the previous section, the operator differential equation H20|l 
with an initial condition H21|l has a solution in the following form 



K(<) - Texp{ / dT [^b-i(T) + b(T)A(r)b-i (r)]} b(a), 
SO we obtain the important identity 

b(i)Texp{ /"*dTA(T)} =Texp{ /* dr [^^b-i(T)+b(T)A(T)b-i (r)]} b(a). 

Ja Ja OT 

(22) 

Now if ^ 

h{t) = Texp{ ( dTB(r)}, 

J a 

then from (|22|l we obtain one of the basic chronological operator identities 
Texp{f dTB(r)}Texp{ / dr A{t)} = 

J a J a 

Texp{ /'*dT[B(T)+Texp{ T B(0} A (t)To exp{- r d^B{0}]}. (23) 

J a J a J a 

If we denote 

C(T)=Texp{/' deB(C)}A(r)Toexp{- / dCB(C)} 

J a 

and solve it with respect to A 

A(t) =Toexp{- r d(B{0}C{T)Texp{ r d^B{0}, 

J a J a 

then it follows immediately from (|23() that 

Texp{/ dT[B(T) + C(t)]} = Texp{ / drB(r)}x 

Texp{ TdrToexpl- T deB(C)} C(t) Texp{ T d^B(e)}}. (24) 

J a 

Identities H23(l and H24|l are generalizations of the well-known Bakcr-Campbell- 
Hausdorff (BCH) and Zassenhaus formulae for "t-dependent" operators in the 
sense that the classical BCH formula merges two exponential operators into a 
single one and Zassenhaus formula splits an exponential operator into a product 
of exponential operators. 

At first sight the identities H23|) and (|24|l represent a circle of some kind. 
To demonstrate the value of these identities we first consider their classical 
expansions and later (in Subsection 5.2 below) we show that in many cases 
of practical important that involve p3() and H24() operator expressions can be 
obtained exactly by a more direct way without using the expansions like ()27f) . 
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3.3 Baker-Campbell-HausdorfF and Zassenhaus formulae 
for chronological operator exponents 

Let us consider the following construction 

K(a)=Texp{/" dr B(r)} A(t) Tq exp{- / drB(r)}. (25) 

J a J a 

If we differentiate K(a) with respect to a we obviously obtain the following 
expression 



5K(a) 
where 



^=Texp{/ dTB(r)}[-B(a)A(i) +A(t)B(a)]Toexp{- / drB(r)}, 

'a Ja Ja 



K{a)\a=t^A{t). 



If we use the commutator notation A(t)B(a) — B(a)A(t) — [A(i),B(a)] we get 
after integration that 

K(a)=A(t)+ pdrTeM T B(0} [A(t), B(t)] Tq exp{- f d^B{0} 

Jt Jt Jt 



or 



Texp{/" drB(r)}A(t)Toexp{- / drB(r)} = 

A(t)- TdrTexpl rdeB(0}[A(t),B(T)]Toexp{- f d^B^O}- (26) 

This is an analog of so-called integral BCH formula. Iterations of leads to 
the following expansion 

Texp{ /" drB(T)}A(t)Toexp{- /" drBir)} ^ 

J a J a 

A{t)~ f dT[A{t),B{r)]+ f dn f dT2[[A(i),B(Ti)],B(T2)]-... 
+ (-1)" f dn f dT2... f dT„ Texp{ f d^B{0}x 
[[...[A(t),B(ri)],B(r2)]...],B(T„)]Toexp{-^*d^B(0}. (27) 
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If we now substitute (|27|l into l|23|l we obtain the more conventional form of 
BCH formula 

Texp{ /" dTB(T)}Texp{ /' dr A{t)} ^ 

J a J a 

Texp{/" dT{A{T)+B{T)~ f dTi[A(r),B(ri)] + 

<j a J a 

f dr^ f dT2 [[A(t),B(ti)],B(t2)] - .. . 

+ (-!)" ( dn I dT2... f dT„Texp{/" deB(0}x 



[[...[A(r),B(ri)],B(T2)]...],B(r„)]Toexp{-^'d^B(C)}}}. 



(28) 



The importanee of the last BCH formula 1)28(1 lies not in the details of the 
formula, but in the fact that there is one, and the fact that it gives the product 
of two operators 

Texp{ /" dTB(T)}Texp{ / dr A{t)} 

J a J a 

as one operator in terms of A and B, brackets of A and B, brackets of brackets, 
etc. 

Analogously, starting from 

K(A) = ToCxp{- / dTB(T)}A(t)Texp{ / dTB(T)} (29) 

<J a <J a 

and differentiating K(A) with respect to A we obtain the foUowing expression 
= To exp{- dT B(r)} [-B(A)A(t) + A(t)B(A)] Tcxp{^^ dr B(t)}, 

where 

K(A)U=a = A(t). 

After integration we get that 
K(A)=A(t)+ f\TToexp{- T B(0} [A(t), B(t)] T exp{ T B(0} 

J a J a J a 

for any A. By setting A = i we arrive to 

Toexp{- / drB(T)}A(t)Tcxp{ / dTB(T)} = 

J a J a 

A{t)+ f drToexpi- f B(0} [A(0, B(t)] T exp{ / d^B{0}- (30) 

J a J a ^ a 
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This is a mirror-like twin brother of the integral BCH formula. We note here that 
is the corollary of ^ if we have in mind the properties of the chronological 
operators ((117)) , lfTT|l and ((T7|) . Iterations of ((ITOI) lead to the following expansion 
(compare with (|27fl ') 

Toexp{- /" drB(T)} A(i)Texp{ /" drBir)} = 

J a J a 

A{t)+ f dT[A(t),B(T)] + / dn / 'dT2[[A(t),B(Tl)],B(T2)] + 

•/a 'J a -J a 

f I dT2... f dT„Toexp{-/" d^B(0}x 

J a J a J a J a 

[[...[A(t),B(ri)],B(r2)]...],B(T„)]Texp{ r d^B{0}. (31) 

a 

If we now substitute H3U|I into 124|) we obtain 
Texp{/" dr[B(r) + C(r)]} = Texp{ /" dr B(t)}T exp{ /" dT[C(r) + 

/ dTiToexp{- rdCB(0}[C(T),B(Ti)]Texp{/' 'd^B(e)}]} (32) 

J a J a J a 

and with the help of the original identity we further obtain 

Texpj/" dT[B(T) + C(T)]} =Texp{ /" dr B(t)} T exp{ /" dTC(T)}x 

J a J a 

Texp{/" dr / dTiToexp{- / C(e)} Tq exp{- / ' B(0} x 

'J a <J a 'J a J a 

[C(t),B(ti)] Texp{ rdCB(0}Texp{ f d^C{i)}}. 

J a J a 

Iterations of (|31() with the help of H23() and H24|l lead to more conventional 
form (but noticeably bulky than (|28|l . so we do not represent it here) of a 
generalized Zassenhaus formula, which nevertheless tells us that the operator 

Texp{/" dT[B(T) + C(T)]} 

J a 

can be expressed as product of operators in terms of A and B, brackets of A 
and B, brackets of brackets, etc. 
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3.4 Relationship between opposite ordered chronological 
exponents and generalized linear operator differential 
equations 

Let us consider now 

Texp{/" dr[A(T) - A(t)]} = 1 = Texp{ / dr A{t)}x 

J a J a 

Texp{- TdrToexpl- T A(e)}A(T)Texp{ T A(e)}}, 

J a J a J a 

SO we deduce at once that 

Toexp{- f drA(r)} = 

J a 

Texp{- / dT To exp{- f d^ A{0}A{t)T cxp{ f d^ A{0}} 



and further 



(33) 



Texp{ f dr A(t)} = 

J a 

Toexp{ fdTToexp{- T A(e)}A(r)Texp{ r d^A{C)}}. (34) 
These identities allow us to solve the following equation 

To exp{ f d^ A(0}A(t)T exp{- f d^ A(0} = B{t) (35) 

J a J a 

with respect to A{t). The solution of (|35|1 is as follows 

A(0 = Toexp{- f d^B{0}B{t)TeM f d^B{0}. (36) 

J a J a 

One can prove this by substitution of into (|^ . Really, in view of and 
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H34|l we have the following chain starting from (|35|) and then substituting 
Toexp{ f d^A{0}A{t)Texp{- fd^AiO} = 

'J a 'J a 

Toexp{/" drToexp{-/" d^B{0}B{T)T exp{ f dCB{0}}x 

J a J a J a 

Toexp{- fd^BiO} Bit) Tcxp{ f d^m)}^ 

Texp{- fdTToexp{- T deB(0}B(r)Texp{ r d^BiO} = 

J a J a J a 

Texp{/" dTB(r)}Toexp{- /" dTB(T)} B(t) x 

J a 

Texp{ /" dTB(T)}Toexp{- / drB(r)} = B(t). 
Let US consider the following operator 

K{t) = To exp{ / dr A(t)} B{t) T exp{ / dr C(r)}. (37) 

J a J a 

Differentiating K(i) with respect to t we have 

dKjt) _ 

dt ^ 

Toexp{^ dTA(T)}[A(i)B(t) + ^?^ + B(t)C(t)]Texp{^ dTC(r)} 
and obtain the following differential expression 

^^^a{t)K{t)+K{t)c{t)+h{t), (38) 



where 



a(i) = Toexp{ /" dr A(t)} A(t) Texp{- /" dTA(T}, 

J a J a 

c(t) =Toexp{- / dr C{t)} C{t) T exp{ f dTC(r)}, 

-/a -/a 

b(t)=Toexp{^*drA(r)} T exp{- £ dr C(t}, 



(39) 
(40) 
(41) 



and moreover 

K{t)\t^a^Bia). (42) 
We can solve (|39|l and H4U|) with respect to A{t) and C(t) respectively 



A{t) = Toexp{- f dTa(r)} a(i) Texp{ /" dTa(T}, 

t/ a J a 



(43) 
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C(t) = -Toexp{ / dr c{t)} c{t) T exp{- f drcir)}. (44) 

J a J a 

From 1)41(1 we have 

^^^^^ =Tcxp{- f drA{r)] h{t) ToGxp{- f drC(r)} 

J a J a 



)} 



dt 

and with H42|l it leads to 

rt PT PT 

B(t)=K{a)+ drTexp{- / A(ti)} b(T) To exp{- / C(ti 

J a J a J a 

and substituting 143|1 and H44() we have, taking into account identities H33|l and 

that 

B{t)=K{a)+ f dr To exp{- f d^ a{m Hr) T exp{- f d^ci^} 

J a J a J a 

If we now express A, B and C in 1)3 7|) via a, b and c then it foUows that 
the formal solution of the operator differential equation (|38|l with the initial 
condition 'K(t)\t=a — K(a), where a, b and c are any linear (not necessarily 
invertible) operators, which do not depend on ^ explicitly, has the following 
form 

K(t) = To exp{ /" dr To exp{- f dn a(Ti)} a(r) Texp{ /" dria(ri}}x 
[K(a)+ f dr To exp{~ f d^ h{r) T exp{~ f d^c{m]x 

Texp{- /" dTToexp{/" c(ti)} c(t) T exp{- /" dric(ri]]. 

J a J a J a 

By using identities (|33|) and H34|) we arrive finally to 
K(t) = Texp{/" dTa(T)}x 

J a 

[K{a)+ f dr To exp{- f d^a(0} b(T) Texp{- / d^c{0}]x 

J a J a J a 

Toexpif drc{r)}, (45) 

J a 

as the formal solution of H38|l . which can be verified directly by substitution of 
^ into (EHl)- 

The equation H38|l is a generalization of equations ^ and 113f) . As we can 
easily see, their solutions © and l(T^ are particular cases of 
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3.5 Differentiation of the chronological exponential with 
respect to a parameter 

In the previous sections we considered some analytical properties of chronologi- 
cal operator exponentials mainly under differentiation on <. It is very important 
from an analytical point of view to consider differentiation of operator exponen- 
tials with respect to some parameter, say a. 
Let us consider the following operator 

^ Texp{^ dTA{T,a)} 

and denote here 

dA(t, a) r d ^ , 

From the identity H26|) we have 

Texp{/" dTA{T,a)} Toexp{- / dTA{T,a)}^ 

Ja OOL Ja 

/ drTexpi A(^, a)} a(T, a) Tq exp{- / d(A{^,a)}, 

J a J T J r 



d_ 

da 



so 



^ Toexp{- / dTA(T,a)} -Toexp{- / dTA(T,a)}[^- 

j drTexp{/" d^A(C,a)} a(r,a) Toexp{-/" d£_ A{£_,a)}]. 

Analogously from identity H3U|) we have 

Toexp{- / drAfr.a)} ^ Texp{ / drA(r,Q;)} = 
Ja oa Ja 

— + dTToexp{- / deA(^,a)} a(T,a) Texp{ / d^A{^,a)}, 

OOL Ja Ja Ja 



which leads to 



3^Texp{/ dTA(T,a)} = Texp{ / dr A{t , a)} [-^ + 
da Ja Ja oa 

f drToexpj- / A(e, a)} a(T, a) T exp{ / d^A{^,a)}]. (46) 

J a J a J a 

3.6 Shift operators 

To transform formal solutions into ordinary expressions we have to express ex- 
plicitly the action of the corresponding linear operators. Every solved equation 
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gives us an example of definite action of a given linear operator. Thus, the list 
of "good" operators is nonempty. Our goal here is to study the properties of a 
certain class of operators that we will refer to as shift operators. 

The Taylor expansion for a sufficiently arbitrary function ^(x) can be ex- 
pressed as 

<i>{x + a)=Y, fci"' d^'J'(^) = exp W^M^) ■ 

k=0 

Here a can be either a constant or a function of arguments which do not include 
X (in the second case ^ is replaced by ^) . If we read this expression from 
the right side, we can obtain the result of the action of an exponential form of 
a shift operator exp{a^} on a function $(a;) (see also (|50ll and (|51|l below): 

exp {a— }$(a;) = $(exp {a—}x) = ^(x + a) 
dx dx 

or 

d d 
exp{a(T) — }cf>(2;) = $(exp {a(T) — }x) = $(a; + q:(t)). (47) 
ox ox 

Let us consider shift operators with a change of variable y — tp{x) . It is 
obvious, that 

exp{a(T) — }4>(y) = exp {a(r)^^}$(V'(a;)) = 
1 d 

exp{a(T);^^ — }$(V(a;)) = $(V'(a;) + ^(t)). 

When ip{x) is fixed function, the set of shift operators with different parame- 
ters a{T) forms an Abelian group. Shift operators with different ^^i^) do not 
commute. 

Simple examples of shift operators are well-known. We give here examples 
which are somewhat exotic 

d 

exp {ln(a)xln(a;) — } $(x) = $(a;") 



and 



exp{^^} exp{ln(2)xln(a;)^} exp{(&-^)^} <S>{x) = ^{x'^ + cx + b) . 

These examples demonstrate that combinations of shift operators may produce 
non-trivial changes of the variable x in the function ^(x). 

The following useful identities for shift operators are known 

exp{a(r).'^^}.^^ exp {-a(r)x"^} = 
exp{a(r)x"^} x^^"" exp {-a(T)x"^} a;"^ = 
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exp{(/3-l)a(T)}x^|; if a = 1; 

2^"^ if " = /5; (48) 

/3(l-e.) Q(l-a) TTT^ c, 

[x 3-° + (1 - a)a{T)x 3-° ] ^ if a / 1 and a ^ (3. 

Alas, in this last type of transformation the "polynomial" nature of the left-hand 
side of the identities in general does not survive in the right-hand side. 



3.7 Chronological operator homomorphism 

As we will see later, the chronological operators with the derivative operator 
A(t) (see definition in Section 2) in the exponential which we will here for short 
denote as ^ 

E = Texp{/ dTA(T)}, 

J a 

play a very important role. 

Let B(i) = b{t) be a function that is the operator of multiplication on the 
b{t). Then, as can be seen from the definition of the derivative operator, the 
commutator [b{t), A(r)] is a function too. 

Now consider the following construction 

K = E6(t)E"\ 

As far as the commutator [&(t), A(r)] is a function, then all repeated commu- 
tators are functions too. So with help of BCH formula H28|) we conclude that 
operator K is a function and 

K = {Eb{t)). (49) 

An analogous conclusion holds for E^^ b{t) E. 

Let us consider further the following obvious chain 

E 6i &2 . . . fo„ = E 6i E"^ E ^2 E"^ . . . E 6„ = 
= (E6i)(E62) ... (E6„). 

With what is observed here and the fact that the chronological operator is 
linear, leads to conclusion that for any function F{bi, 62, . . . , &«), which can 
be expanded in power series with respect to foi, 62, we can obtain the 

following nice property: 

E F{bi, 62, ... , 6„) = F((E6i), (E62), . . . , (E6„)). (50) 

Analogously 

E-i F(6i, 62, ... , 6„) - F{{E-^ 61), (E~i 62), . . . , (E-^ 6„)). (51) 
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3.8 Some additional useful identities 

By taking inverses of both sides of (|23|l wc find the similar identity for the 
product of inverse ordered operators 

Toexp{/" dTB(T)}Toexp{ / dr A{t)} ^ 

J a J a 

Toexp{ rdT[B(r) + Tcxp{- T B(0}A(t)To exp{ f di^^S,)}]]. 

•J a J a J a 

Furthermore, the analog of (|24|) is 
Toexp{/" dT[B(T)+C(r)]} = 

J a 

Toexp{ fdrToeM T deB(0}C(T)Texp{- TdeBlOllx 

J a J a J a 

Toexp{/ dTB(T)}. (52) 

J a 

If in H22|l we substitute 

b(t) = To CXp{ /' dTB(T)}, 
^ a 

then from (|22|l we obtain the following operator identity 
Toexp{/" dTB(T)}Texp{ /" drA(T)} = 

Texpl/ drToexpif d^B{0} [A{t) +B{t)]T exp{~ f d^Bi^]}, 

J a J a J a 

and by inverting it and making the change of operators A — s- — B and B —A 
we obtain another form 

Toexp{/" dTB(T)}Texp{ / dr A(r)} = 

Toexp{- TdrToexpl- T A(e)} [A(r) + B(r)]Texp{ T A(0}]}. 
Since from (|24|) we have that 
Texp{/" dT[B(T) + C(r)]} = Texp{ / dTB(T)}x 

Texp{ f dT To cM- rdCB(0}C(T)Texp{ r d^Bi^} = 

J a J a J a 

Texp{ /" dTC(r)}x 

^ a 

Texpif dr To exp{- f C(0} B(t) Texp{ / dSC(0}}. 
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then we come to a third form 



Toexp{- / dTC(r)}Texp{ / drBir)} ^ 

J a J a 

Texp{ f dT To eM~ T C(0} B(t) T exp{ r dCCiO}}^ 

J a J a J a 

Toexpl-/ dT To exp{- f d^B{0} C{t)T exp{ f d^B{0}}. 

J a J a J a 

Let US now consider the invertible operator B which does not depend on t. 
Then from H22|l we have 

BTexp{/" dr A{t)} = T exj>{ f dTBA{T)B-^}B. 

J a J a 

or ^ ^ 

Texp{/ dTA(T)}B = BTexp{ / dTB"^A(T)B}. 

./a 

From the expression 

K(i) =BTexp{ /" dTA(r)B} 

^ a 

we obtain by differentiation that 

9K(t) _ g ^^^^ g Texp{ / dr A(t)B} = B A{t) K(i), 



(53) 



dt 



K(t)|t=, -B. 



Solving the above operator differential equation we find that for any linear (not 
necessarily invertible) "i-independent" operator B 

BTexp{/" dTA(T)B} = Texp{ /" dTBA(T)}B. 

J a J a 

Analogously, 

BToexp{/" dTA(T)B} = ToCxp{ /" drBA(T)}B. 

J a J a 

3.9 Example 1. Linear first-order PDE 

The general linear first-order partial differential equation for u{t, p) 

du{t,p) , ^ , . /, ^ 1^ ^\ ^ f ^ du{t,p) : f ^ du(t,p) 

— Q-^ — ^ ni.P)+h{t,P)u{t,p)+fi{t,p) —^— + - ■ ■ + J,,{t,p) , (54) 
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where we abbreviate p = xi, ...,a;„, with initial condition 

is fitted outright to be solved by the operator method. Its formal solution follows 
from ifTHl) : 

u{t, p) = Texp{^ dT [/o(t, p) + /i(r, + • • . + J^{r, p)-^]} v{p)^ 



£ dr T exp{^* di P) + h{Lp)-^^ 



We indicate that the above general solution can be expressed in terms of 
some particular solutions of homogeneous first-order PDE: 

If we introduce the following notation 



(t, p) = T exp{y dr [/i (r, p) — + ■•• + /„ (r, p) —]} : 



(z = l,. 



we can consider it as the fundamental set of particular solutions of H55() with 
initial conditions Ci{t, p)\t=a — Xi. Since the operator in the exponential is the 
derivative, then taking into account (|5()|l . the solution of homogeneous first-order 
PDE 1)55(1 may be expressed via Q as 

u{t,p) = v{(i{t,p), ...,Cn{t,p)). 

The last expression is the general solution of the homogeneous first-order PDE 
H55(l if v{p} is an arbitrary function. 
If we further denote 

/■* 9 d 
Z,{t,T,p) = Texp{J dT[fi{T,p) — ^ h/„(T,p)— , (i = l,...,n), 

we can find by using ((Tn|l that 
or 

Cl{t, p) = Ci{T, Zl{t,T,p),..., Zn{t, T,p)), 

SO Zi{t,T,p) are solutions of the following system of algebraic equations 

Q{T,Zi,...,Zn) = Ci{t,P), {i = l,...,n). (56) 
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Let us take up in passing a nice property of considered type of functions 
which we will need later. Denoting 

/■* a d 

h{t,p) = To exp{~ dr [/i(r,/9) — +■ • •+/„(r,p)— ]}a;j , (i = l,...,n) 

(57) 

then from and H51|) we have 

Toexp{- / dT[hiT,p)— + --- + UT,p)—]}Q{t,p) = 

Qit,biit,p),...,bnit,p))=x, (58) 

and analogously 

bi{tXi{t,P),---,Cn{t,p))=Xi, (l = 1, . . . ,7l). (59) 

That is (i and bi are bundled by algebraic systems (|58ll and (|59ll . 
Moreover, it is obvious that 

= -To cxp{- £ dr [/i(r, p) A + . . . + /„(^, /^(^, , 

and from (|51|) it follows that 

^^^ = -/,(i,6i(t,p),...,6„(t,p)), (z = l,...,n). (60) 

The expressions (|57ll are solutions of the system of n first-order non-linear ODEs 
(EOI). 

As by virtue of and lO 

dr [Mr,p) + hir,P)-g^^ + ' ' ' + /^(t,^)— ]} 



dXn 

T exp{ /* dr [A (r, p) ^ + ■ • • + /„ (r, } x 

exp{£ dr (To exp{^' [/i($, + . . . + /„(^, f^A_]}f^^r, p))} = 

exp{ / dr fo{T,Zi{t,T,p),...,Znit,T,p))}x 

J a 

Texp{ r dr [h{r,p)^ + ■■■ + fn{r,p)j^]} . 

Ja UX\ OXji 

If we now return to the general linear first-order PDE (|54|l . we can rewrite its 
operator solution via (and Zi which are expressed through by system (|56|) ') 
into the following form 

u{t,p) = u(Ci(i,p), ...,C«(i,p))exp{ / dr fo{T,Zi{t,T,p),...,Zn{t,T,p))}+ 

J a 

/ dT(j){T,Zi(t,T,p},. . . ,Z„(i,T,p))exp{ / d^fo{^,Zi(t,S„p), . . . , Z„(t, ^, p))}. 

t/ a J T 
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3.10 Example 2. Linear parabolic differential equation 

The linear parabolic differential equation 

— g-^ — = Jo[t,P)u[t,p) + /i(t,p) \ ^ fn(t,p) —g^ — ' 

where p = xi, ...,a;„, with an initial condition 

u{t,p)\t=o = v{p), 
has the following formal solution 



/■* 52 
P) = T exp{y dT[/o(r,/3) + /i(T,/9)^ 



2 52 

•• + •/'"(^'^a^2"]}"(^■ 



2 
1 

It is interesting to note that it is generally accepted that the parabolic equa- 
tion is a second order PDE, but from the angle of the operator method con- 
sidered here this equation is first order if the initial condition is specified for t 
variable. 

In the simplest case when all fi are constants fi — ki then it is easy to see 
with (|SSl that 

/■* 92 g2 
Texp{ / dr [ko + ki—^ H h fcn^-j]} exp{i{aiXi H h cr„a;„)} = 

Jo ^^n 

exp{i(CTia;i H h cr„x„)} exp{t [fco - kial fcncr^]} , 

so we can obtain well-known non-operator solutions of such a parabolic equation 
with help of a Fourier transformation of the initial condition. 



4 Non-linear first-order ODE's 

4.1 The operator solution for first-order non-linear ODE 

An ordinary non-linear first-order differential equation 

du{t) 



dt 



^fit,u{t)) (61) 



on account of its non-linearity is not immediately suited for application of the 
above considered operator method. Nevertheless there are many possibilities to 
convert the problem to a linear one and in Subsection 3.9 we have obtained 
by way of the operator method a solution for such an equation. Here we consider 
this problem in more detail. 

Most linearization procedures are concerned with the introduction of spaces 
of larger dimensions. Here we demonstrate some of them which are almost 
generally applicable. 
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If we introduce a new function 

S-C^, w) = e'^"(*\ (62) 

from one can easily derive the following first-order (with respect to t) linear 
differential equation for S{t,u!) equivalent to 

^^L!^^u;f{t,^)S{t,cu) , S{a,c.) = e-<'^\ (63) 
from which the formal solution follows immediately from in the form 

Sit,,.) = T exp{£dra;/(r, ^)}e'^"(°). (64) 
So the solution of equation l|6HI is: 

.(t) = ^Uo. (65) 

Direct substitution of the solution (|65|) reduces the equation to an iden- 
tity. 

The theorem for the uniqueness of a solution of a differential equation shows 
that if we have an exact solution of a problem in several forms, all the forms 
can be transformed to each other. So the main problem that remains in this 
approach is how to transform the operator solution to that form which can be 
considered most useful. 

4.2 Solution for first-order ODE in form with derivative 
operator 

The main result of this subsection is another form Hfi7(l of the formal solution for 
non- linear first-order ODE. We can prove the solution correctness by its direct 
substitution into ODE Ht)l|) . but it is desirable to outline the way which leads 
to such form. 

Let us start from an operator solution of the equation H63|l in the form H64|l 
but let us rewrite it (not only) for convenience in the following notation 

S{t,c,uj)=T exp{J^ dTUjf{T,-^)}e'^', 

where c = u{t)\t=a- With the help of identity H53|l we find 

S(t, c, uj) = e"^" T exp{ f druj e-"7(T, T^)e'^n ■ 1. 

Ja OUJ 

CHB-expansion gives 
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Since it is easily proven by induction that 

,Lu\ ^ m 



.m — 1 



then from expansion of f{x) into power series and backwards summation it 
follows that 

therefore we can conclude with taking into account l|()t)|) and the property of the 
shift operator H47|) that 

/(,,c+A)=„p{A|)/(r,c)»p{-A|,^ 

Then 

5(<,c,a;) = e-=exp{^^}T exp{£ dr - exp{-|-|-} ■ 1 



e"^exp{A^}T exp{^*dr [a;/(r,c) 



9c ' 9cj 9c 



By expanding the chronological exponential of the derivative operator /(r, c)^ 
with the help of identity (Pljl and properties (|5n|l we obtain that 

S{t, c, Lo) = e^^^ exp{^ ^ }T exp{- j'^ dr /(r, c) |; } x 

exp{- / dT5(T,c)}exp{w / dTG{T,c)}, 

J a <y a 



where 



and 



g{T,c) ^ To exp{J\u{^,c)-^J 



d,df{T,c) 



dc 



G(r,c) = To eM[dU{^,c)^^fiT,c). 
If we differentiate the last expression with respect to uj wc find that 



dS{t, c,cj) 
duo 



c S{t, c, uj) + 



d d /■* d 

^^^Pi^TrlT exp{- / dT/(T,c) — }; 
aw oc oc 

exp{— / dr (7(r, c)} expjcj / dTG{T,c)} I dTG{T,c). 

J a J a J a 
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If in the second item we recover the initial operator form 
dS{t,c,Lj) 



duj 



c S{t, c, UJ 



) + Texp{^ dTLof{T,-^)}e'-- dTG(r,c), 



and note that the chronological operator commutes with any function which 
does not depend on cj, wc arrive at 

dS{t,C,Uj) , , f\ M r./, X 



du) 

and inasmuch as S{t,c,u})\i^=Q = 1, then 



u(t,c) = — |^^o = c+/ dTG[T,c) 



hence it follows that 

u{t, c) — C' 

or finally 



u{t, c) = To exp{^ dT f{T, c) ^} c. (67) 

The expression H67|l is the general solution of the equation (|61|l if we will consider 
c as an arbitrary constant. 

Above derivation is not rigorous, of course, so it is very important to verify 
our conclusion by direct substitution of the obtained solution into equation (|61(l . 
After differentiation of (|67|l and using property H5()|l we can be sure that H67|l is 
really the formal solution of the equation H61|l . 

So we have obtained a solution for a non-linear first-order equation in a more 
convenient form with the derivative operator. 

Let us finish this Subsection with a brief remark about the well-known in- 
terconnection between first-order ODEs and linear first-order PDEs. 

As we have seen in Example 1 above (see Subsection 3.9) the function 

C(i,c) =Texp{-£dT/(T,c)^}c (68) 

satisfies the equation 

^ + /(i,c)^=0. 
ot oc 

If we now act on both sides of (|68|) by the operator 

To exp{^ dr/(r,c)^}. 
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we get 



C{t,u{t,c}) = c 



(69) 



and 

u{tX{t,c))=c 

- the well-known equations expressing the relationship between solutions of 
first-order ODEs and the corresponding linear first-order PDEs. 

In some practical cases it is easier to obtain explicit expression of (^{t, c) 
for given ODE, for example, by an integrating factor method. (I69II then is the 
implicit form of the solution for the given ODE. Which form of solution of types 
(|67|) or (|69|) is more analyzable depends on the object under investigation for a 
specific problem. 



4.3 Example 3. Bernoulli ODE 

The formal solution of a Bernoulli ODE 

— = a{t)u" + b{t)u 
in the form of 167|l is as follows 

u{t, c) = To exp{ / dr [a(r)c" + 6(t)c] —} c. 

Jo oc 

Expanding the above chronological exponential with the help of identity H52|) 

u{t, c) = To exp{ j dT a(r) exp{ J d^ KC)c ^} c" ^ exp{- d^ ^(^)^ 

/■* d 
exp{ / drbMc—Jc 
Jo 9c 

and with H48|l gives us the solution with only shift operators 

u(t,c)—exp{ dr aM expUa — 1) / d(, b{^)} c°' —} exp{ drbMc—jc. 

Jo Jo oc Jo dc 

Executing shift operations we arrive to a classical form of the solution for 
Bernoulh ODE (a ^ 1) 

u(t,c) = [c(i~"' + (l-a) / dra(r) exp{(a-l) / d^b{0}]^ exp{ drbir)}. 



The technique being used here is applicable not only for Bernoulli ODEs, but 
it is successful under concatenation of some circumstances, when the original 
chronological exponential is decomposed into the chain of shift operators on its 
own without solving any auxiliary differential equations. 
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4.4 The "integral-free" form of solution for non-linear first- 
order ODE 



Let us consider the following operator chain 

d , EH 



exp{{t~a)[f{s,c) — + —]} = 



To exp{ J dTe«^''^)*}/(s,c)^e«''-^);t^>} exp{(t-a)^} 

/■* d d 

To exp{ J dTf{s + T-a,c)—} exp{{t-a) — }. 



So 



/■* d d d 

To exp{ j dTf{s + T- a, c) — } c = exp{(t - a) [f{s, ^)'q^ + g^]}^ 



and if we denote 



d d 

U{t,c,s) = exp{(t - a) [/(s,c) — + ^]}c 



we obtain from the preceding expression that the general solution of ODE H61|l 
is 

U(t,c) = U{t,C,s)\s=a 

or 

d d 

u{t,c) = [exp{(i - a) [/(s,c) — + —]] c\\s=a ■ 

The last form of the operator solution has some interesting features in the 
sense that it does not contain any integration and even ordering operator T, 
which may be useful for calculating approximate expressions of the solution 
uit, c). 

The function U {t, c, s) obviously satisfies the following PDE 
dU{t,c,s) , dU{t,c,s) dU{t,c,s) 

-^r- - ^) dT- ^ ° ' ^' ^)"*='^ = 

4.5 The solution with an arbitrary function 

Let us now consider f{t, c) = g{t, c) + h{t, c) so we can find the following form 
of the general solution of equation (|()T|l 

d d 
c) = To exp{y (iT/(T,c) —}c = To exp{y dr [g{T, c) + h{T, c)] —} c ^ 

To exp{^ dT[To exp{^ d^h{^,c) ^}5(t,c)] x 

exp{-/^de[Toexp{^^.CMCc)|}^]}|} 



X 



/■* d 
To exp{ J dr hij, c) — } c 
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or 



u{t, c) = To exp{^ dT [To exp{^ h{^, c) ^} [/(r, c) - h(T, c)]] x 

exp{- y d^ [To exp{ dC h{C, c) ^} ■ 

/■* 9 
To exp{ J dr h{T, c) ^} c, 

where we can consider h{t, c) as an arbitrary differentiable function. 
Supposing it is known that 

/■* d 

z{t,c) = To exp{ J dTh{T,c)—}c, 

then with the help of property (|50|l we obtain 

u{t,c) = To eM J dri^^^£^^^Jz{t,c). (70) 

As far as h{t, c) is an arbitrary function, the preceding expression is vahd for 
any differentiable function z(t, c). The particular solution of equation with 
initial condition u{t, c)\t^a = c is obtained when z(t, c)\t=a — c. 
Rewriting ITUJ as 

u{t, c) = z(t, To exp{ dr i^^±j->L^^ ^} c) (71) 

dc 

by wisely choosing the form of the function z(t,c) we may reduce the initial 
problem to 

u{t, c) = To eM J dr ——^^—}c, 

dc 

which may be a simpler one by losing some troublesome singularities of the 
initial problem. 

5 The systems of non-linear first order ODEs 

5.1 The formal solution of the system of non-linear first 
order ODEs 

To find the formal solution of the system of non-linear first order ODEs 

^^ = /.(t,ui(t),...,u„(t)), (z = l,...,n) (72) 
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we, of course, could introduce an auxiliary function like and pass through 
a chain of unwieldy expressions but now we are ready to assert that the general 
solution of the system (|72|l has the following operator form 

=Toexp{ / dT[/i(T,c)^ + --- + /„(T,c)^]}c,, (73) 

Ja OCi dCn 

where we abbreviate c = ci, c„ and Ci are a set of arbitrary constants. 
Differentiating (|73|l with respect to t and with (|51|) we have 

5^ = To exp{ / dr [/i(r, c)^ + ■■■ + /„(r, c)^]} /.(t, c) = 

at Ja OCi dCn 

= fi{t,Ul, . . ■,Un). 

If we denote here that 

C^(t,c)=Texp{- / dT[/i(T,c)— + --- + /„(t,c)— ]}0(c), 

Ja OCi dCn 

it is easy to see that for any function satisfies a PDE similar to H55|) : 

So if we know the n fundamental solutions of this equation, namely 

I — 

do. 



0(i,cO = Texp{- / dT[fi{T,c}^ + --- + U{T,c)^]}c,, 



then 

C0(i,5) = 0(Ci(i,5),...,C«(i,^) 
and we can find Ui as the solution of the algebraic system 

Ci{t,ui, . . . ,M„) = Ci , [i = 1, . . . ,n). 

Analogous to the first order ODE, we can rewrite the solution 173|) in an 
"integral-free" form as follows 

d 
u,{t,c) = [exp{(t - a) [fi{s,c)- h • • • + /„(s,c)-- h ^]} c,]\s=a ■ 

Oci OCn Os 

5.2 The direct calculation of BCH type expressions when 
the involved operators are derivatives 

As we have seen, the operators similar to ()29|l when operators A{t) and B(t) 
are derivatives similar to play a key role in our approach. And one of the 
principal points here is the calculation of the following expressions 

K,{t,x) = 

Toexp{-^ dTj2hj{T,x)-^}^Texp{J^ dr J^hjir, x)-^} , (74) 
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where x — xi, . . . , Xm, hj{t, x) are ordinary functions, and we suppose that for 
a given hj{t, x) we are able to calculate the fundamental set of of functions 



/■* " d 
Zj(t,x) = Toexp{- / dTy^hj{T,x)- — } x^ 

Ja t"^ dXj 



We have calculated above such expressions for m = 1, but if m > 1 we need to 
follow by a different way. 

Analogous to if^ . we find that 

dK,it,x) 



dt 

r , v"^ , / 5 , dhjit.x) d , 

Toexp{-y^ ..X:.,(r,f)-}(X:^^-)x 

Texp{/ dT^/i,(r,f) — }. 

The use of H51|l leads to the following system of linear operator ODEs 

^ f2 g,, (t, x)K, it,x), = 1, . . . , m) , (75) 

where _ 

d dhj{t, x) 



f d 

gij{t, x) = To exp{- dr ^ h.j{T, x) — } 



dxi dxi 
with initial conditions 

K,{t,x)\t=a = -^. (76) 
We can write down the formal solution of the system (|75|) - (|76|1 as 

/t rn m q 

i=ij=i ' 

Since the system (|75|) is linear, then its solutions depend on initial conditions 
linearly, i.e. 

™ d 
K, (<, f ) ^ Pik (t, x) — , 

fe=i 



where 

Pik 



dr '^'^gij{T,x)cj jr-} Ci]\ c^=i.c,^^^o ■ 
1=1 j=i "^^ 



This conclusion is a direct consequence of the fact that the set of operators 
12]Li hj{t, x)-£p forms a Lie algebra. 
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So in cases under consideration we can exactly calculate the operators 'Ki{t, x) 
ifTHl without resorting to BCH expansions. 

The usage of such calculations in principle allows us to factor the chrono- 
logical exponential, e.g. in (|73|l . into a product of relatively simple factors (with 
a sequence order assigned in advance) 

Toexp{/ dT^/j(r,c) — } = [] To exp{/ dr g,{T,c) — } . 
Unfortunately this way leads to complicated non-linear PDEs for functions 



6 The non-linear nth order ODE 

For the non-linear nth order ODE 

dTuit) duit) d^-VO 

we can apply the results of the previous Section as long as the nth order ODE 
can be represented by a system of n first order ODEs: 



dt 

duj (t) _ 
dt ~ 

dUn-l(t) 

dt 



Ut+i{t); 

= f{t,u{t),Ui{t), . . .,Un-l{t)). 



So the formal solution of ODE H77|) can be expressed in the following form (here 

C — Ci , . . . , Cji ) 

d d d 

dr [fir,(^-Q^ + ■■■ + + ■■■+ C2^]} ci , (78) 

and it is obvious that (here m = l,...,n — 1) 

d'"u(t,c) ^ r f\ ,r, ^ d d d 

= To exp{ I dr [/(r, c)^ + • • • + C+i^ + • • • + c^^]} c,„+i . 

If we consider ci, c„ as arbitrary constants, the expression H78|) is the general 
solution of the equation H77|l . 

As we have indicated in Example 1 (see Subsection 3.9), the functions 

/■* a d d 

G(t,c) = Texp{- / dT[/(T,c)- \ hCi+i— H hC2-^]}ci, 
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obey the differential equation 



-K7 + f[s,c)- h--- + c,+i— H hc2-^ = [Q\t^a = ct) (79) 

Ot OCn OCi OCi 

and 

' ■ ■ ■ ' ) = ^' ' (* = 1> • ■ • > ri). 

So if we know some, say k < n, independent solutions of H79|) . then we can 
eliminate k unknowns from the system 

Qit.u,—,..., )^c^, {t = l,...,k). 

and as a result reduce the order of initial problem from n to {n — k) . 
In the "integral-free" form 

d d d d 

u{t,c) = [exp{(i - a) [f{s,c)— H ^^'-+'^0^. ^ ^^^9^ ^ Ss^^*^^^ ' 

the auxiliary function 

C/(t, c» = exp{(t - a) [/(., c) A + . . . + + . . . + + ^]} 

satisfies the following PDE 

du ^ac/ ac/ ac/ du 

and 

Since here we have a larger number of degrees of freedom than in the case 
n = 1, we can obtain many equivalent forms of the formal solutions like H71|l 
using the way described in the previous section. 



7 Helmholtz equation 

Here we consider the essential features of the operator method for linear par- 
tial differential equations on example of formal solution of Helmholtz equation 
under different formulations of the boundary conditions. The solutions of the 
Helmholtz equation represent the (spatial part of) solutions of the wave equa- 
tion. 

Let us consider the Helmholtz equation for an inhomogeneous medium, which 
has the following form 

d u{x,y,z) ^ _^ ^^^^ z)]u{x, y, z) + q{x, y, z) (80) 
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with boundary conditions on x — a 

'U'{x,y,z)\x=a^ a{y,z) , = /3(y, , (81) 

where A2 = ^ + e{x,y,z) is a function, which takes into account the 
wave speed dependence on the space point, and q{x, y, z) is a wave source. 
If we introduce the function 

Q( u \ u I \ I du{x,y,z) 
b(x,y,z;h,p) = hu{x,y,z) +p — , 

where h and p are auxihary real parameters, then equation 18U|I can be trans- 
formed into the following form 

^'^'"^'fe^'^'^'* " ~ ['^2 + £ix,y,z)]p-^]S{x,y,z;h,p) +pq{x,y,z) 

with boundary condition on a; = a 

S{x,y,z;h,p)\^=a = ha{y,z) +p/3{y,z) . 

Hence the solution of the Helmholtz equation (|80|l under boundary conditions 
(El) is (a; > a) 



d r d d 

u{x,y,z) = — Texp{ j dT[h— - [A2 + e{T,y, z)]p —]} [ha{y, z) + p f3{y, z)]- 



^ ^ dr T exp{^ ^ " +sity,z)]p-^]}p q{T, y, z) . (82) 

It is easy to see that u{x, y, z) in H82I) does not depend on auxiliary parameters 
h and p. 

There are boundary value problems for Helmholtz equation (which are more 
profound from a physical point of view) when one puts certain requirements on 
the solution behavior at infinity (at r = {x^ + 2/^ + z'^Y^'^ — > 00). Since here 
the co-ordinate r is selected by boundary conditions it is expedient to solve 
this problem in spherical co-ordinates with the initial supposition that u and 
its first derivative on r at r = a are known. At a ^ 00, u and its derivative 
in a medium with Ime > have to diminish rapidly enough (there are not 
wave sources at infinity), hence the first item of the expression of type 
goes to zero. Therefore for u to satisfy the radiation conditions at infinity in a 
boundless inhomogeneous medium we can find the following operator expression 
(/to £ > 0) 

8 r°° r h r) F B c 

-W = -^ / 'iCToexpiy^ ^C[^^-pA + ^M^1]p^]}K^.(^f)), 

(83) 

where vector r = {x,y,z) and A — -\- + is the ordinary Laplace 
operator. 
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The Green function for a boundless inhomogeneous medium follows from 
H83() at q{r) = 5{r — tq). By using the well-known method of images one can 
find from H83|) the Green function say for half-space a; > and to get a solution 
of a Helmholtz equation (|80|l one has to fulfil the radiation conditions at infinity 
and a boundary condition on a plane for example. 



8 The system of linear first-order PDEs 

Let us consider the systems of linear PDEs, which are first-order with respect 
to X and y, using as an example the following two PDEs for one function 

u{x,y,p) = u: 



^ = Aix,y)u, 



du 

— = B{x,y)u, 
dy 



(84) 



where p is a set of parameters say zi, z„, A{x,y) = A{x,y, p) and B(a;, y) 

d_ 

dy 



B{x,y,p) are linear operators, which do not depend on ^ and ^ explicitly. 



We can solve the first of them by obtaining 

= Texp{ / dTA{T,y)}(j){y,p), 

J a 

where (/)(y, p) is a yet unknown function. Substituting now this solution into 
the second equation of the system we have 

^Texp{^ dTA{T,y)}^{y,p)^B{x,y)T exp{J^ dr A{t , y)} (jy{y , p) 

and with the help of l|46(l we obtain the differential expression: 

= [To exp{- £ dT A(r, y)] B(x, y) T [ dr A(t, y)]- 

r dr Toexp{- f A(g, y)} T exp{ f A(^, y)] ] 4>{y. p) , (85) 

J a J a J a 

which have the form of differential equation for an unknown (j){y, p) if and only 
if the right-hand side of (|85|l does not depend on that is when its derivative 
with respect to x is equal to zero, which leads to the well-known consistency 
condition 

\\( \ Tx( M , ^H^^y) dB{x,y) 

[A(x,y),B{x,y)] + — — — = 0. (86) 
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We can solve with the help of H45|l this operator equation with respect to 
B(a;, y), that is rewrite the consistency condition as 

B(a;,y) = Texp{ / dTA{T,y)}x 

J a 

[B{a, y) + 1^ dr Tq exp{- d^ A{t v)} T exp{- d^ A(e, y)}] x 

Toexp{- / dTA{T,y)}, (87) 

J a 

where B(a,?/) — Ti{x,y)\x=a is an arbitrary linear operator (which docs not 
depend on ^ and ^ explicitly). 

The formal solution of the equation 18511 is 

py px px 

</.(y,p) = Texp{ / dCToexp{- / dr A(r, ()} B(a;, C) Texp{ / drA(r,C)}- 

J b J a J a 

dr To exp{- £ d^ A{t ()} ^^^^ T exp{^" d^ A(^, C)}} c{p) , (88) 
where c{p) is an arbitrary function. Substituting now (|87f) into H88f) we receive 
</.(2/,p)=Texp{ r dC'B{a,C)}c{p) 

Jb 

and finally the desired solution of considered system H84() if (|86() holds is 

u(x,y,p) =Texp{ / dr A(t, y, p)}Texp{ T dC B(a, C, p)} c(p) , (89) 

Ja Jb 

where c{p) has an obvious definition as c{p) = u{x,y,p)\ 

Since here A{x,y,p) and 'B{x,y,p) are operators, then H89(l can represent 
solutions of non-trivial systems of PDEs. 



9 Conclusions 

We have presented some ways in solving DEs by the chronological operator 
method. Besides linear first-order DEs and systems of such DEs, we have ob- 
tained operator solutions for linear and non- linear ODEs of arbitrary order. 

It is easy to note that the obtained solutions contain differential operators 
with respect to arbitrary constants, which represent initial conditions of the 
problem. For more complicated problems, e.g. for non-linear PDEs, the formal 
solutions will contain variational differential operators. Some examples we have 
touched on can be found in 9 and 10 . 

In conclusion, we believe that we succeeded in demonstrating the fact that 
operator forms of DE solutions can be handled analytically no worse than ordi- 
nary functions. In some cases its transformation properties are more comfortable 
than, e.g., for some special functions. 
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